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■ Abstract 
^ ■ 

ly-^ I Numerical integration formulas in n-dimensional Euclidean space of degree three are 

■ discussed. For the integrals with permutation symmetry we present a method to construct 
its third-degree integration formulas with 2n real points. We present a decomposition method 

' . and only need to deal with n one-dimensional moment problems independently. 

■ Keywords: Numerical integration; Degree three; Cubature formulae; Decomposi- 
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1 Introduction 

^ ■ Let = R[ xi, . . . , Xri\ be the space of polynomials in n real variables and ^ be a square 

QQ , positive linear functional defined on H" such as those given by -Sf (/) = f{x)W{x)dx, where 

OO ' is a nonnegative weight function with finite moments of all order. Let Yi^ be the space of 

in 



o 



X 



polynomials of degree at most d. Here we discuss numerical integration formulas of the form 

j^{f)^Y.^kf{u^'^), (1) 



where are constants and u^^'^ are points in the spaces. The formulas are called degree of d if 
they are exact for integrations of any polynomials of x of degree at most d but not d+\. 

In this paper, we only deal with the construction of third-degree cubature formula. It 
■ looks like a simple problem, however it remains to be solved. The well known result is due to 

Stroud [ll-[3] • He presented a method to construct numerical integration formulas of degree 3 for 
centrally symmetrical region and recently Xiu [5] also considered the similar numerical formulas 
for integrals as 



^(f) = ... I w{xi)...,w{Xn)f{xi,X2...,Xn)dxi...dx. 



n- 



In 0], Xiu assumed that every single integral is symmetrical, which means his result naturally 
belongs to centrally symmetrical case. Recently the authors Q extended Stroud's results and 



*The project is supported by NNSF of China(Nos. 61033012,11171052,11271060,61272371) and also supported 
by "the Fundamental Research Funds for the Central Universities" . 
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presented formulas of degree 3 of 2n points or 2n + 1 points for integrals as 



^{f)=j ... j Wi{xi) . . .Wn{Xn)f{xi,. . . ,Xn) dxi. . .dXn, 

Wi{xi)>^, Xj G [aj,6i], i = l,...,n. 



Besides, many scholars employed the invariant theory method to deal with symmetrical case and 
we can refer to [THSll and the references therein. As far as we know, 2n is the minimum of the 



integration points except some special regions(see [Ul]), and for centrally symmetrical region 
Mysovskikh had shown this point. However, for the general integration case, it remains 
unknown how to construct the formulas of degree 3 with 2n points. In the two-dimensional 



case, third degree integration formulas with 4 real points was given in Ij, [ij] for any regions. 



But it is difficult to extend it to higher dimension. For other related work, we can refer to 



15l4l7l| and the reference therein. 

This paper will extend the results in for the integrals of product regions to those with 
permutation symmetry. First we present a condition which is satisfied by the integral. And then 
we prove that under this condition, the construction problem of cubature formulae with degree 
three can be transformed into two smaller sub-cubature problems. Finally, for the construction 
of cubature rules of the integrals with permutation symmetry can be decomposed into n one- 
dimensional moment problems. 

This paper is organized as follows. The construction of cubature formulas of degree 3 are 
presented in section 2. And section 3 will present two examples to illustrate the construction 
process. Finally, section 4 will make a conclusion. 



2 The construction of third-degree formulas 

Assume that ^ has the following property: 

There exist n linearly independent polynomial li{xi, . . . , x„), i = 1, 2, . . . , n 
(P) such that all liln{i = 1, ... ,n — 1) are the orthogonal polynomials of degree 

two with respect to S£ . 

Let 

T : li{xi, . . . ,Xn) — >ti, z = l,2, ...,n (2) 

be a linear transformation and ^ be transformed into »Sf' . Then by the assumption all titn{i = 
1, 2 . . . , n — 1) are the orthogonal polynomial of degree two with respect to Here we do not 
require that all titn{i = 1, 2 . . . , n — 1) can constitute a basis of orthogonal polynomials of degree 
2 with respect to . We can also assume that the third-degree formula of .if' has the following 
form 

v'-^^ = (^^l,l,1'l,2,'yi,3 • • . ,'Ul,n-l,0) UJl 
7;(2) = (U2,l, i;2,2, ^2,3 • • • , t^2,n--l, 0) UJ2 

(3) 

t)W = {VN,1,VN,2,VN,3 ■ ■ ■ ,VN,n~l,0) ^^N 
^(N+j) = (0, 0, . . . , 0, VN+j,n) i^N+j, j = 1, 2. 

To enforce polynomial exactness of degree 3, it suffices to require ([3]) to be exact for 

1) ^1 ) ^2, ■ ■ ■ ,tn, tjtj , titjtj; i, j, k — 1, 2, . . . , 71. 
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Then we have 



OOi + UJ2 + . . . + UJN+2 = ^'(1) 
f UJlVi^i + UJ2V2,i + . . . + iONVN,i = -^'{k), i = 1,2, ... ,11 - I 
\ ^N+lVN+l,n + ^^N+2VN+2,n = {tn) 

( UJiVi^iVij + UJ2V2,iV2J + . . . + UJNVN,iVN,j = ^'{Utj), i,j = 1,2, 
I '^TV+l^'ir+i,™ + UJN+2V%^2,n = i^l) 

f iOiVi^iVi^jVi^k + • • • + (^NVN,iVNjVN,k = ^'{Utjik), k = 1,2, 

and the equation (|4ap can be rewritten as 

UJl + U>2 + ■ ■ ■ + UJN = £.1, 
+ '^N+2 = S,2 

6 + e2 = ^'(i). 

Hence we can rewrite the equations (HD as 

UJi + U}2 + ■ ■ ■ + I^N = (l, 

UJiVi^i+UJ2V2,i + ■■■+^NVN,i= ^'{ti), i = 1,2, . . . ,n - 1 
UJlVi^iVij + UJ2V2,iV2,j + • • • + UJNVN,iVN,j = ^'{Utj), i,j = 1, 2, 
^ UJlVi^iVijVi^k + • • • + UJNVN,iVN,jVN,k = ^'{Utjtk), k = 1,2, 

UJN+l + = ^^2 

i^N+lVN+l,n + <^N+2VN+2,n = (tn) 
WiV+lW^+l,„ + UJN+2V%^2,n = ^'(*n) 
. WiV+l<+l,„ + UJN+2V%+2,n = -^'(4) 

6 + e2 = i?'(i). 

Once is determined by ([Sc]) . then ()5ap and ()5bp become one n — 1 dimensional and the 
other one-dimensional moment problems respectively. If these two lower dimensional moment 
problems can be solved, then we can get a cubature formula of degree 3 with respect to the 
original integration problem. Generally speaking, the one-dimensional moment problem can be 
easily solved, but it is difficult to be solved for the n — 1 dimensional moment problem. However, 
if the n—1 dimensional problem can be divided into one n — 2 dimensional moment problem 
and the other one-dimensional moment problem and further the n — 2 dimensional moment 
problem can continue this process, then the original integration problem can be turned into n 
one-dimensional moment problems. 

In what follows, we will prove that if ^ is a integral functional with permutation symmetry, 
then the construction problem of third-degree cubature formulae can be turned into n one- 
dimensional moment problems. In fact, we usually encounter this kind of integral functional, 
for example, the integration over the simplex, the square, the ball or the positive sector of the 
ball, that is {{xi,X2, ■ ■ ■ , > 0,i = 1,2, . . . ,n; xf + x"^ + . . . , < 1}. 

We first prove that the integral functional with permutation symmetry must meet the prop- 
erty (P). Thus the original cubature problem can be turned into two sub-cubature problems. 

Theorem 1. Let n >2 and 

^2 (-^1 5 • • • 5 '^n) — "^i ^ — li^,... , n 1 

Inixi, . . . ,Xn) = Xi+ X2 + . . . + Xn + Cn 

where 

^{x\ + (n — 3)xfx2 — (n — 2)xiX2X^) 

" .^{x\ — X1X2) 



,n — 1 
,n — 1 



(4a) 
(4b) 

(4c) 

(4d) 



1 (5a) 
. . . , n — 1 

. . . , n — 1 

(5b) 
(5c) 
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If ^ is permutation symmetrical, then liln (i = 1, 2, . . . , n — 1) is orthogonal to the polynomials 
of degree < 1 . 

Proof. Take liln as an example. We first prove .^{x\ — X1X2) 7^ to confirm the existence of 
In- In fact, by the symmetry and the positivity, 

^(xf - X1X2) = ^^{xf - X1X2 - X1X2 + xl) = ^^{{xi - ^2)^) > 0. 
Let us exam the orthogonality of liln- Assume n > 3. By the symmetry, we have 

.^{hln) = ^{x1- xl + (Xi - Xn){x2 + • • • + Xn-1 + C„)) = 0, 

-^{Xihln) = ^{x\xi - x\xi + (Xl - Xn)(x2Xi + . . . + X^-lXj + C„Xi)) =0, 2 < i < U - 1 

and for i = l(the same to i = n) 

=5f (xi/l/„) = ^{xf - X^Xl + (Xi - X„)(X2X1 + . . . + Xn-lXl + C„Xl)) 

= ^(xf + (n - 3)x^X2 - (n - 2)xiX2X3 + c„(xf - X1X2)) = 0. 
It is easy to verify that the result holds when n = 2. This completes the proof. □ 

Let and be two linear functionals defined on 11""^ and 11^, whose moments are 
determined by 

^i{t,t,tk) = ^{hijik), ^liutj) = ^{kij), ^i{ti) = ^{k), ^i(i) = er> (6) 

+ = ^(1) and l<i,j,k<n-l 

respectively. Thus the construction problem of third-degree formulas with respect to ^ is turned 
into two smaller problems, one of which is the construction with respect to =Sfi and the other of 
which is the construction with respect to »Sf^. It is easy to compute 

Jfi(t„) = if = n ■ ^(xi) + c„^(l), 

^^{tn^) = n^(x\) + n(n - l)^(xiX2) + 2nc„ • ££(xx) + c^^(l) 
Jfi(t„3) = n^{x\) + 6Q-^(x?X2) + 6Qif(xiX2X3) 

+3c„(n^(x?) + n(n - l)^(xiX2)) + 3nc^^(xi) + 4^(1) 
= nJf (xf ) + 3n(n - l)^(x?X2) + n(n - l)(n - 2)^(xiX2X3) 
+3cn{n^{xl) + n(n - l)^(xiX2)) + 3n4^(xi) + 4^(1). 

Next we will show this decomposition process can continue. 

Let us consider the cubature formula with respect to Jfi. Obviously, _Sfi is also permutation 
symmetrical, which allows us to employ theorem [1] continuously. Define 

4^\'ti, ■ ■ ■ ,tn-l) = t'i — tn-1, i = 1, 2, . . . , n — 2, 
In-lih, . . . ,tn) = ti + t2 + . . . + tn-1 + Cn-1, 

where 

_ ^i(tf + (n - 4)tjt2 - (n - 3)tit2t3) 
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then «f^^2i,n = 1,2,... ,n — 2 arc the orthogonal polynomials of degree two with respect to 
=^1. Noticing ti = li{xi,X2, ■ ■ ■ , Xn), we have 

(^1; • • • 1 ^n—l) ~ ti in— 1 ~ 3?^— 1, i = 1,2, ... ,71 2, 

4-1 (^1' ■ ■ ■ ,tn) = ti + t2 + . . . + tn-1 + Cn-1 

= XI + X2 + ... + Xn-1 - (n - l)Xn + C^,-!. 

Again let ^2 and .=Sf^ be two linear functionals defined on n""^ and H^, whose moments are 
determined by 

.^2iT^TJTk)=mPlfl^^), ^2(x.r,)=^i(/f)/f ), ^2 (t,) = =^1 (if ^ =5f2(l) = 



2 



^?+d'^=^^i^ and l<z,i,A;<n-2 
respectively. It is easy to compute 

^2(r„_i) = ifi(/«i) =^tf2i) + c„-i(^i(l) -^(1)), 

^^Tn-l^) = = ^1 ([(*1 + ^2 + . . . + tn-l) + C„-l]') 

= ^(tf2,)2) + 4_l(^l(l)-^(l)), 
^2(r„_i3) = + - ^(1)) . 

Assume that jSf^ is a linear functional defined on I!"-*^ for every k{0 < k < n) and satisfies 
property (P). Then a cubature problem of degree 3 with respect to Jffe can be divided into 
two smaller cubature problems — one with respect to ^k+i and the other with respect to 
Moreover ^k+i also satisfies the property (P) and then this process can continue and will end 
when k = n — 1. Finally, an n-dimensional cubature problem can be transformed into n one- 
dimensional cubature problems. 

Theorem 2. Let 

Inlk = Xi+X2 + ...+ Xn-k - {u - k)Xn-k+l + Cn-k 

and be a linear functional defined on U.^ according to the above process, then the correspond- 
ing moments are 



^>^+\r) = ^{ll^l,) + cn-k{^k{l)-ni)) 

= C„_fc^fc(l), 

= (n - k){n -k + - X1X2) + ^.^^^(l), 

^fc+l(^3) ^ if((/('^:),)3) + ct,(^,(l)-^(l)) 

= (n - k){n -k + l)(n -k + 2)if (-x? + 3x^X2 - 2x1X2X3) + 



where 



^{xf - X1X2) 
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Proof. It remains to prove Eq.(l8|). It follows from theorem [T] that 

_ ^fc(ti + (n - 3 - k)tlt2 - (n - 2 - k)tit2h) 

According to the definition of , we have 

■^k{t\ + (n - 3 - k)t\t2 - (n - 2 - k)tit2tz) 
= ^{tl + (n - 3 - k)tlt2 - (n - 2 - k)tit2h) 

= ^{{Xi - Xn-kf + (n - 3 - k){xi - Xn-kf{x2 - Xn-k) 

-{n - 2 - k){xi - Xn-k){x2 - Xn~k){x2. - Xn-k)^ 

= ^( — 3x^X2 + 2x1X22:3 



^k{tl-tit2) = ^{{Xi - Xn-kf - {Xl - Xn^k){x2 - Xn-k)) 
= J^'{xl - X1X2), 

where the permutation symmetry is used. This completes the proof. □ 

Remark 1. In fact J^n-i is also a one- dimensional integration functional. The corresponding 
moment can be calculated by 



^n-l{l)=it"\ 

^„_i(r) = if(xi -X2) =0, 

^n-l{r^) = ^{{Xl - X2)2) = 2if (X? - X1X2) 
^„„l(T3)=if((xi-X2f)=0. 



(9) 



For convenience, in what follows let = J^n-i- 
Suppose that 



if '=(5) « ^iMti,k), A: = 1, 2, . . . , n (10) 

i=l 

is exact for any (7 € Ilg. And let v*''^ = (1^^ "Wj- • • • , vf^k) be the solution of 

Xi + X2 + . . . + X„ + Cn = 

Xji—1 x^i — 

for /c = 1, (11) 

X2- Xn = 
^ Xi- Xn = 

Xi + X2 + ■ ■ ■ + Xn + Cn = 

+ X2 + . . . + Xn-1 - {n- l)Xn + C„_i = 

for A; = n, (12) 

2^1 + X2 - 2X3 + C2 = 
^ Xl X2 — ti,n 
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and 



Xi + X2 + ■ ■ ■ + Xn + Cn = 

X1+X2 + ■■■+ Xn-1 - (n - l)Xn + Cn-1 = 



X1+X2 + ■■■+ Xn-k+1 - {n-k + l)Xn-k+2 + Cn-fc+1 = U^k 
•^n—k •^n—k+1 



for A; = 2, 3, . . . , n — 1, 



^ Xi - Xn-k+l = 

Hence the final cubature formula can be written as 

n rih 



k=l 1=1 



which is exact for any / E Ilg. It is clear that the solution of Eq. ()lip is 



{r]i,r]i,. . . ,r]i), r]i 

and the solution of Eq.(|13p is 

Cn-l — Cn — 52hti2 



n 



Xr, 



n 



_ Cn-2 — Cn-l _ 

n — 1 

Xn-k+^ = Xn~k+i 

Cn—k+2 A;+3 ^i,k ti,k 
Xn-k+2 — Xn-k+Z H ; — TTi — Xn-k+Z 



Xi = X2 



n-k + 2 

Xn-k+l = Xn-k+2 + 



ii,k C-n—k+l 



n-k + 2 



n — k + 1 

where 52, fc = 1 if = 2 and 52,k = if /c 7^ 2, and the solution of Eq. (fT2]l is 

Cn— 1 Cn 



Xn — Xn—l — • • • — X^ 

C-n—l ~l~ ti n 
X2 = X^ 



n 



(13) 



71 ^71 — 1 

Xi=X2+ ti^n = X3 + — 

Collecting the above discussion, we have 

Theorem 3. Assume that ^ is •permutation symmetrical. Then there must exist a cubature 
formula 



mi 



■^U) ~ '^'^Wi,kf{ai,k, ■ ■ .,ai^k,Pi,k,l,l,-- . ,7) + • • • '"i.i) •= ^(f) 

_ . ^ I ^ NX 



k=2 i=l 



n-k+1 



k-2 



i=l 



which is exact for every polynomial of degree < 3. In the formula, a 's and j3 's can he computed 



7 



by 



7 



Cn-l 



n 



Pi,k + 



n — k + 1 



2<k<n-l; 
2<k <n-l: 



n. 



(3i,n = 7 



ti,n + C2 



a,; 



where ti^j. 's and Wi^k 's are the nodes and weights of the quadrature formula (|1U|) with respect to 
The proof is a direct result of the computation and is omitted. 

Remark 2. For the quadrature problem of the one- dimensional moment, it is well known 
that the number of the nodes nj- = 2 in the general case. Hence the total number of the nodes of 
the cubature formula with respect to J£ is 2n generally and 2n is usually the minimum among 
all the cubature formula of degree 3 except one case of the integration over the n- dimensional 
simplex fTjj /. For more knowledge of the problem of the one- dimensional moment, we can refer 
to the appendix of Q/. 

Remark 3. For convenience, we present the relations of ^'s as follows 



^(1) 



+ 

?2 



+ 

?2 



(n-2) 



+ 

Jn-2) 
?2 



According to the previous discussion, it is clear that 
.(k-i) 



^fc(l) = for l<k<n-l and ^"(1) = = ^1 



(n-2) 



3 Numerical Examples 

• Firstly take the integration over the n-dimensional simplex as an example. Define 



■^(/) = / f{xi,X2, x„)dxid2;2 . . . dxn 



where 



Tn = {ixi,X2, . . . ,Xn)\xi -\- X2 + . . . + Xn < 1, Xj > 0, i = 1,2, 

It is well known that 



\Xi -^2 ■ ■ ■ -^n 

Then by a simple computation, we have 

2 



ai!Q2! . . . Otn\ 



(n + Qi + Q2 + . . . + a„)! ' 



n + 3 
1 



i — 2,3,. ..,72 1, Cfi 



n + 2 
"n + 3' 



n + 3 



,n}. 
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Here if we take 



and 



and 



S2 



t 



i+l 



(n-2) 



t 



n • n! n • n! 

n, then the moments of ^^^^{1 < k < n — 2) are 



, i = 0, 1, . . . , n — 2, 



n ■ n\ 



Cn-k-^k(X) 



(k) ^2^ 



2(n-Eti4) ^ 
(n + 3)n • n! ' 

+ c2_fc(^fc(l)-^(l)) 



n 



A:)(n-fc + l) , 4(n-Etitfc) 



+ 



^n~kJ J ' -n-k\-^y-J -v-/y- (n + 2)! ' n(n + 3)2 • n! ' 

(n - A;)(n - + l)(n - + 2)Jf (-x? + 8x1X2 - 2x1X2X3) + ^.^^^(l) 
-2(n - A:)(n -k + l){n-k + 2) 8{n - J^Li ^k) 



(n + 3)! 



n{n + 3)^ • n! 



n • n\ 



I' 



(n+l)!(n + 3)' 

+ 5n + 8 
(n + 3)!(n + 3)' 

2(n + 2)(n + 4) 
~(n + 3)2(n + 3)! 



^"(1) 
^"(r) 



n • n! 



I' 



0, 



(n + 2)!' 



0. 



By taking different values for ^s, we can get different cubature formulae. For n = 3 and n = 4, 
if we take all tj = 1, then we can get formulas as showed in Tables [T] and [2j 



Xi 



X2 



X3 



weight 



0.34240723692377 
0.14125289379518 
0.41353088165296 
0.12380973765487 
0.60719461208592 
0.05947205458075 



0.34240723692377 0. 

0.14125289379518 0. 

0.41353088165296 0. 

0.12380973765487 0. 

0.05947205458075 0. 

0.60719461208592 0. 



34240723692377 
14125289379518 
00627157002742 
58571385802358 
16666666666667 
16666666666667 



0.01469064053612 
0.04086491501944 
0.01887111233337 
0.03668444322218 
0.02777777777778 
0.02777777777778 



Table 1: Nodes and weights for T3 

In Tabled! the first point is outside of the region. To avoid it, we can take 

_ 93 _ 378 _ 108 
^ ~ 85' ~ 391' ~ US' 
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Xl 

0.27145760185760 
0.12024746726682 
0.30652570925957 
0.11154151763119 
0.37131176827505 
0.09266869570542 
0.54391317546145 
0.0275 1539596712 

Wl 



£i 

0.33333333333333 
0.14285714285714 
0.41666666666667 
0.12037037037037 
0.61593041596355 
0.05073625070311 



X2 

0.27145760185760 
0.12024746726682 
0.30652570925957 
0.11154151763119 
0.37131176827505 
0.09266869570542 
0.02751539596712 
0.54391317546145 



X2 

0.33333333333333 
0.14285714285714 
0.41666666666667 
0.12037037037037 
0.05073625070311 
0.61593041596355 



X3 

0.27145760185760 
0.12024746726682 
0.30652570925957 
0.11154151763119 
-0.02833782226438 
0.52894832287488 
0.14285714285714 
0.14285714285714 



£3 

0.33333333333333 
0.14285714285714 
0.00000000000000 
0.59259259259259 
0.16666666666667 
0.16666666666667 



X4 

0.27145760185760 
0.12024746726682 
-0.06243427063585 
0.52251830424930 
0.14285714285714 
0.14285714285714 
0.14285714285714 
0.14285714285714 



weight 
0.01875000000000 
0.04203431372549 
0.01875000000000 
0.03495843989770 
0.02608695652174 
0.02608695652174 



= 0.00254167472911, W2 = 0.00787499193755, W3 = 0.00294495824332, 
= 0.00747170842335, W5 = 0.00365639117145, wq = 0.00676027549522, 
= 0.00520833333333, ws = 0.00520833333333. 



Table 2: Nodes and weights for T4 



Table 3: Nodes and weights for T3 



and the corresponding cubature formula is listed in Table [3j In the formula, the first and third 
nodes are on the boundary of the region T3. 
If we take 



94 

85' 



t2 = 1, h 



76 

85' 



then all the nodes are inside the region, see Table HI 



Xl 



X2 



X3 



weight 



0.33237874197689 
0.14303319621635 
0.41247250927755 
0.12085734331926 
0.61593041596355 
0.04371016618787 



0.33237874197689 
0.14303319621635 
0.41247250927755 
0.12085734331926 
0.04371016618787 
0.61593041596355 



0.33237874197689 
0.14303319621635 
0.00838831477823 
0.59161864669481 
0.16666666666667 
0.16666666666667 



0.01927056497746 
0.04216734351927 
0.02043165185637 
0.03512390369919 
0.02516339869281 
0.02516339869281 



Table 4: Nodes and weights for 

In Table El all the weights are positive. However, we find there exist three points outside of 
the region. If we take 

_ 104 _ 3577 _ 9947 _ 49 
~ 75"' ~ 2775' ~ 8880' ~ 60' 



and add one more point with weight —49/80, then the corresponding formula is showed in Table 

El 

In Table O there are 5 points on the boundary. To make all the nodes inside the region, we 
can take 



7 187 

5' 145' 160283 



179522 5 



and add one more node with weight —618391/961698 and the corresponding formula is shown 
in Table El 
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Xl 

0.25000000000000 
0.12500000000000 
0.28571428571429 
0.11224489795918 
0.35714285714286 
0.07792207792208 
0.00000000000000 
0.42857142857143 
0.28571428571429 



X2 

0.25000000000000 
0.12500000000000 
0.28571428571429 
0.11224489795918 
0.35714285714286 
0.07792207792208 
0.42857142857143 
0.00000000000000 
0.28571428571429 



X3 

0.25000000000000 
0.12500000000000 
0.28571428571429 
0.11224489795918 
0.00000000000000 
0.55844155844156 
0.14285714285714 
0.14285714285714 
0.14285714285714 





0.25000000000000 
0.12500000000000 
0.00000000000000 
0.52040816326531 
0.14285714285714 
0.14285714285714 
0.14285714285714 
0.14285714285714 
0.14285714285714 



wi = 0.00555555555556, W2 = 0.00888888888889, W3 = 0.00600490196078, 
W4 = 0.00742227521640, W5 = 0.00633074935401, we = 0.00533755957993, 
W7 = 0.00425347222222, ws = 0.00425347222222, lUg = -49/80. 



Table 5: Nodes and weights for T4 



Xl 



X2 



X3 



X4 



0.24955035825246 
0.12511271452382 
0.28567845223177 
0.11210697374487 
0.35714280321315 
0.07749399446444 
0.56855699818890 
0.00287157323967 
0.28571428571429 



0.24955035825246 
0.12511271452382 
0.28567845223177 
0.11210697374487 
0.35714280321315 
0.07749399446444 
0.00287157323967 
0.56855699818890 
0.28571428571429 



0.24955035825246 
0.12511271452382 
0.28567845223177 
0.11210697374487 
0.00000010785942 
0.55929772535683 
0.14285714285714 
0.14285714285714 
0.14285714285714 



0.24955035825246 
0.12511271452382 
0.00010750044754 
0.52082193590823 
0.14285714285714 
0.14285714285714 
0.14285714285714 
0.14285714285714 
0.14285714285714 



wi = 0.00566710734383, W2 = 0.0089162259895080, ^3 = 

W4 = 0.00739793083678, W5 = 0.0063627387083462, wa = 
W7 = 0.00434027777778, wg = 0.0043402777777778, wg = 



0.006035977209200, 

0.005780572945942, 
-0.643019950129875. 



Table 6: Nodes and weights for T4 



The integration problem on the n-simplex was studied very extensively. According to the col- 
lection of R. Cools in the website(http://nines.cs. kuleuven.be/research/ecf/ecf. html), the mini- 
mum number of nodes in the third-degree formulas is n + 2, in which there is a negative weight. 
Except the (n + 2)-point formula, the minimum number is 8 and 10 for n = 3 and n = 4, 
respectively. If we only consider the formulas with positive weight, the minimum of the points 
is 8 and 11 for n = 3 and n = 4 respectively. Therefore our formula for n = 3 and n = 4 have 
the fewest numbers among all the formulas with positive weights. 

• Secondly take the integration over the positive sector of a ball as an example. Define 



^{f)= / f{xi,X2,...,Xn)dxidX2...dXn 

" Sn 



where Sn = {(xi, X2, . . . , Xn)\xf + + . . . + < 1, xi > 0, . . . , > 0}. It is easy to get that 
where Uo denotes the number of the odd number among ai, a2, • • • , and mU denotes the 
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double factorial of m and m!! = 1 if m < 0. Then by a simple computation, we have 



Ci = - 



(n + 2)!! /7rJV]-[f] 4- 



(n + 3)!! V2 



IT 



7r-2' 



2,3, ...,n-l, 
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(n + 2)!! /7rN[=fi]-[f] (n - l)7r - (2n - 4) 
(n + 3)!!V2y ' tt - 2 

(n + 2)!! /7rx[^]-[t] 



(n + 3)!! V2 



Here if we take 



S2 



H+1 
n ■ n\\ 



TT\ [§] Jn-2) 
' SI 



n ■ n!! 



2, 



0,1,2, 



,n-2. 



then the moments of ^'^"'"^(1 < k < n — 2) are 
^^=+^(1) 



.77.1! *^2^ ' 



n • n 



(n + 2)-(n-Etiii) 



n- (n + 3)!! 

(n-/c)(n-A: + l) /7r\[^] _ -,\ 
(n + 2)!! V2J ■ V2 ~ J 



(n + 2)! 

(«-Ef=ii.)[(« + 2)!!]' /4-7rx2 .j,.2[^]-[^] 



n-nW- [(n + 3)!!]2 



TT-2 



(n- A;)(n- fe + l)(n- fc + 2) /TrV"^! /vr 



(n + 3)!! 



(?^-Ef=i*i)[(" + 2)!!]2 M-ttnS /7rN3[^]-2[t] 



n-n!!-[(n + 3)!!]3 



/4 — TTN-l /Try 

V7r-2J ' V2y 



and 

^i(l) = 

^\t) = 

^\t^) = 



n ■ nil 



2, 



2(n7r - 3n + TT - 4) /vry^] 



(n + 3)!! ■ (7r-2) 

n(f + n-l) /7r\[^] 2n-(n + 2)!! 
(n + 2)!! V2J 



+ 



(n + 1)!! • (n + 3)!! V2 
7r\2[V]-[f] /(n-l)7r-(2n-4)\2 



TT-2 



(n + 2)! 

n + 2 
[(n + 3)!!]2 

3n^-n /TTxlVl n(n- l)(n-2) /7r\[^] 
(n + 3)!! ■ V 2 J (n + 3)!! V 2 J 

3n(f + n-l) ^^y^l (n - l)7r - (2n - 4) 

.2, 



7r-2 



7r^2[ii^]-[f] (n - l)7r - (2n - 4) 

)^ 



(n + 3)!! 
3n 



TT-2 



(n + 1)!! L(n + 3)!! 



(n + 2)!! 



TT 



3[ 



n — 3 1 or ^ 



■{n-l)TT- (2n-4)\2 
. 7r-2 / 



1 



(n + 2)!! i3 /7r\3[^]-2[t] / (n - l)7r - (2n - 4) - 
(n + 3)!!J ' V2 J ' V tt - 2 . 
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and 



If take all ti = 1 for n 
E If we take h = 0.8, t2 = 



^"(1) = 
^"(r) = 

= 3 and n 
1.31,^3 = 



n ■ n\\ 

0, 

2 

(n + 2)!! 
0. 



[f] 



vr 

1 

2 



= 4, then we can get formulae as showed in Tables [7] and 
1.11 and ^4 = 0.78 for n = 4, then we can get a formula 



with all the nodes inside the region, see Table [9l 



Xl 



X2 



weight 



0.53887049476004 
0.18341741723402 
0.57520979290336 
0.20283315000517 
0.76016315955181 
0.09981758853698 



0.53887049476004 
0.18341741723402 
0.57520979290336 
0.20283315000517 
0.09981758853698 
0.76016315955181 



0.53887049476004 
0.18341741723402 
0.02206116228206 
0.76681444807844 
0.31250000000000 
0.31250000000000 



0.07852747507104 
0.09600545012840 
0.06975676243570 
0.10477616276373 
0.08726646259972 
0.08726646259972 



Table 7: Nodes and weights for S-^ 



Xl 

0.47721483105875 
0.17126237887529 
0.48420041705925 
0.20526869095372 
0.56004995494835 
0.16531879543241 
0.74847573599445 
0.05241038692400 



X2 

0.47721483105875 
0.17126237887529 
0.48420041705925 
0.20526869095372 
0.56004995494835 
0.16531879543241 
0.05241038692400 
0.74847573599445 



£3 

0.47721483105875 
0.17126237887529 
0.48420041705925 
0.20526869095372 

-0.02818760532449 
0.76127471370739 
0.29102618165375 
0.29102618165375 



X4 

0.47721483105875 
0.17126237887529 
-0.06966276495181 
0.76713241336478 
0.29102618165375 
0.29102618165375 
0.29102618165375 
0.29102618165375 



wi = 0.03771636146294, W2 = 0.03938992292057, = 0.02874740384082, 
WA = 0.04835888054269, = 0.03167997303102, wq = 0.04542631135249, 
w-j = 0.03855314219176, 0.03855314219176. 



Table 8: Nodes and weights for ^4 



4 Conclusion 

In this paper, we present a method to construct third-degree formulae for integrals with 
permutation symmetry. Our method is a decomposition method, which is easy to compute. At 
the end of the paper, we present some numerical results, which seem to be new. Compared 
with the existing method, we focus on the case of permutation symmetry, which seldom is 
considered. The numerical results show that the number of the points attain of close to the 
minimum. Besides, in most cases, the weights in our formulas are all positive or at most one 
negative weight. 
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